Magnetic Percolation and the Phase Diagram of the Disordered RKKY Model 
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We consider ferromagnetism in spatially randomly located magnetic moments, as in a diluted 
magnetic semiconductor, coupled via the carrier-mediated indirect exchange RKKY interaction. We 
obtain via Monte Carlo the magnetic phase diagram as a function of the impurity moment density rii 
and the relative carrier concentration n c /ni. As evidenced by the diverging correlation length and 
magnetic susceptibility, the boundary between ferromagnetic (FM) and non-ferromagnetic (NF) 
phases constitutes a line of zero temperature critical points which can be viewed as a magnetic 
percolation transition. In the dilute limit, we find that bulk ferromagnetism vanishes for n c /ni > .1. 
We also incorporate the local antiferromagnetic direct superexchange interaction between nearest 
neighbor impurities, and examine the impact of a damping factor in the RKKY range function. 
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There has been substantial recent interest in the old 
problem of long-range ferromagnetic ordering in localized 
"impurity" magnetic moments induced by indirect ex- 
change ("RKKY") interaction mediated through (effec- 
tively) "free" carriers (either electrons or holes) . The re- 
cent interest arises from the context of ferromagnetic or- 
dering in diluted magnetic semiconductors (DMS), such 
as Gai-zMn-jAs with x « 0.1, where the Mn dopants 
act both as the impurity magnetic moments and as 
acceptors producing the carriers (which happen to be 
holes for GaMnAs) mediating the RKKY coupling. The 
standard model for DMS ferromagnetism has been the 
carrier-mediated RKKY interaction, and understanding 
the RKKY magnetic phase diagram, therefore, takes on 
particular significance. This is important in view of the 
diluted and the random nature of the spatial distribu- 
tion of the impurity magnetic moments which could lead 
to substantial frustration in the magnetic interaction be- 
tween the impurity moments due to oscillatory nature of 
the RKKY coupling. The latter yield substantial anti- 
ferromagnetic (AF) couplings which have the potential 
to disrupt ferromagnetism. Our goal in this paper is to 
numerically obtain the RKKY magnetic phase diagram 
in a disordered DMS system via direct Monte Carlo sim- 
ulations. 

The complex interplay of the the long-ranged oscilla- 
tory behavior of the RKKY interaction and strong dis- 
order makes simple theoretical statements difficult, and 
it is not obvious a priori whether in the dilute limit the 
ferromagnetic state is supported for any choice of sys- 
tem parameters even at T — 0. Mean field treatments 
such as the Curie Weiss continuum mean field theory 
(MFT) are problematic, because they fail to take into 
account the discrete crystal lattice and thermal fluctu- 
ations and assume a ferromagnetic ground state with- 
out providing any means of assessing the validity of this 
assumption. In this Letter, we rigorously take into ac- 
count positional disorder and we obtain the true ground 
state spin configurations, finding that the RKKY model 
does support a ferromagnetic phase, albeit only for a 



limited parameter range. In addition, we demonstrate 
via direct Monte Carlo that the transition from the non- 
ferromagnetic (NF) to the ferromagnetic (FM) phase is 
marked by the percolation of magnetic clusters which 
grow in size as AF couplings are reduced and ultimately 
coalesce to yield long-range ferromagnetic order; this con- 
stitutes zero temperature percolation critical behavior. 
We neglect all quantum fluctuations, but our interest be- 
ing the interplay of disorder and magnetic interaction, 
our results should in general be valid with respect to the 
existence (or not) of the FM phase. 

Although the specific calculations described in this 
work are motivated by possible carrier-mediated RKKY 
ferromagnetism in DMS systems, our results also ap- 
ply more broadly to a variety of disordered magnets 
with competing interactions (e.g. the Edwards- Anderson 
model which we have also examined and found similar be- 
havior) where the NF to FM transition at (T=0) occurs 
via magnetic percolation. Note that we do not attempt 
to characterize the NF phase (which may be a simple 
paramagnet or a subtle glassy phase) except to empha- 
size that it does not have long-range FM ordering. 

Our model physical system is Gai-^Mn^As (by far 
the most studied DMS system) where Xi ~ 0.01 — 0.1, 
the concentration of Mn dopants, is in the dilute limit. 
We assume that Mn impurities only occupy Ga sites in 
the zinc-blende GaAs (fee) lattice with a lattice con- 
stant a. The large spin (S = 5/2) of the Mn moments 
permits a classical treatment of the spins, which we re- 
gard as Heisenberg spins governed by the Hamiltonian 
W = Eij A r «)Si-Si+E(tj> J^Si-Sj, where in the first 
term is the separation between moments i and j and 
J(r) is the RKKY range function given in three dimen- 
sions by J(r) = Joe~''/ i r~ 4 [sin(2/cF7') — 2kpr cos(2kpr)]; 
kp = (§7r 2 n c ) 1 / 3 is the Fermi wave number, and n c is the 
hole density (a closely related quantity is x c = n c /4, the 
number of holes per lattice site). Jq(> 0) is related to the 
local Zener coupling J p d between the Mn local moments 
and the hole spins. In particular, we have Jo oc mJ^ d 
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with m being the hole effective mass and the S 2 — (5/2) 2 
factor being absorbed into J . With DMS materials be- 
ing at best poor metals, we also introduce in the RKKY 
range function a damping factor e~ r l l where the damping 
scale I is related to the carrier mean free path or perhaps a 
carrier localization length arising from Anderson localiza- 
tion of the free carriers; however, we initially set I = oo in 
order to study the long-ranged RKKY coupling. For the 
latter, an important subtle question is whether the full 
RKKY model, with its long-ranged oscillatory behavior, 
supports a ferromagnetic ground state for strongly disor- 
dered systems. In the second term of H, J AF is the local 
antiferromagnetic superexchange coupling which is rele- 
vant only for neighboring impurities on the fee lattice. 
We note that very large values of Jo (which tend to lo- 
calize the hole carriers and hamper indirect exchange [1] ) 
are also deleterious to ferromagnetism and lower T c , but 
this effect is not examined in this paper. We also ignore 
all band structure effects, which should be adequate for 
qualitative purposes. Our goal here is to obtain explicitly 
the ideal RKKY phase diagram, rather than calculate re- 
sults for a specific material or compare with experiments. 

Salient length scales include the mean spacing between 
impurities, l s = n\ ^ 3 and the scale of the oscillations 
in the RKKY range function, fcp 1 . Hence kpl s deter- 
mines the relative importance of ferromagnetic and an- 
tiferromagnetic couplings; for kpl s ~ 1 RKKY oscilla- 
tions tend to disrupt ferromagnetic order yielding a NF 
ground state, while one expects ferromagnetic interac- 
tions to dominate for fcp/ s <C 1. Since kpl s oc (nc/ni) 1 / 3 , 
and due to its experimental relevance, n c /rii is a useful 
parameter of merit for gauging the importance of AF in- 
teractions. While we concentrate on ferromagnetic order- 
ing for n c /rii <C 1 regime, other work has examined the 
n c /rii 3> 1 limit, deep in the NF phase [2,3]. We will see 
that, although a stable ferromagnetic phase is supported 
by the pure RKKY model, it occurs only for a relatively 
narrow n c /n,i domain; in particular, in the dilute limit, 
we find long-range ferromagnetic order at T = only 
for n c /rii < .1. Previous Monte Carlo calculations [4-6] 
have examined on a qualitative basis the impact of com- 
peting interactions on the ferromagnetic state. We show 
explicitly that the transition of the ferromagnetic to the 
non-ferromagnetic phase at T = is a disorder driven 
critical transition involving magnetic percolation. 

In the presence of strong disorder, n c /rii has a role very 
similar to temperature (k^T / Jq); just as thermal fluctu- 
ations at higher temperatures disrupt ferromagnetic or- 
der by flipping spins, a larger n c /rii is associated with 
strong AF interactions which prevent many pairs of spins 
from aligning. For large enough n c /ni, only spins in close 
proximity are ferromagnetically correlated. As n c /rii is 
decreased, these small correlated groups of spins increase 
in size. Ultimately, the magnetic clusters span the en- 
tire system and magnetic percolation occurs, signaling 



the appearance of long-range ferromagnetic order. This 
is the extended free-carrier analog of the bound magnetic 
polaron percolation ferromagnetic transition recently dis- 
cussed [7] in the context of strongly localized DMS mate- 
rials. To identify magnetic percolation and thereby locate 
the NF/FM phase boundary, we determine the typical 
size of the magnetic clusters using a standard technique 
to calculate the ferromagnetic correlation length £ within 
Monte Carlo [8]. Magnetic percolation and concomitant 
long-range ferromagnetism occurs when £ becomes com- 
parable to the system size L. 

In our Monte Carlo calculations, we average over 
at least 500 disorder realizations. We obtain ground 
state configurations via Monte Carlo simulated anneal- 
ing; thermal fluctuations are provided by the Heat Bath 
technique [9]. To exploit finite size scaling, we exam- 
ine the behavior of the normalized correlation length £/ L 
(where L is the linear dimension of the system) as a func- 
tion of system size. In the NF phase, £/L diminishes 
as L is increased, while for ferromagnetic order, £/L in- 
creases with increasing L. One seeks the FM/NF phase 
boundary where £/ L is constant in L for the critical value 
of n c /rii. To minimize finite size effects, we examine a 
range of system sizes where the mean number of spins 
(N) contained in the system is at least on the order of 
1000. To calculate the Curie Temperature T c , we also 
examine £/L, but we vary the temperature T instead of 
n c /ni. Again, in obtaining T c , we seek the critical £/L 
curve. 

In Fig. 1, we display calculated characteristics which 
would be accessible in experiment; these results are ob- 
tained for a large system with (N) ~ 1000. Panel (a) 
of Fig. 1 is a graph of the ferromagnetic order param- 
eter m = [(|S|)] (square brackets indicated disorder av- 
eraging), the magnetization obtained for the undamped 
(I = oo) case. For small n c /m, one sees what ap- 
pears to be a plateau where spins are essentially per- 
fectly collinear; a reduction in the polarization begins for 
larger n c /rii with the magnetization becoming strongly 
attenuated for n c /rii ~ 0.1. It is tempting to regard the 
"plateau" and noncollinear regimes as indicative of qual- 
itatively distinct phases, but this notion can be seen to 
be illusory if one considers that for any spin, there is al- 
ways a finite probability of having a void large enough 
that the distance from the spin to its nearest neighbor is 
greater than the distance to the first zero of the RKKY 
oscillation. Since in this situation the moment would not 
interact ferromagnetically with its nearest neighbor, it 
is reasonable to assume that there is always some non- 
collinearity for any finite value of n c /rii. 
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FIG. 1. T = magnetization and susceptibility, and T c 
plots versus n c /n; for Xi = 0.1; results are shown for the 
undamped RKKY with no AF superexchange coupling. Pan- 
els (a) and (b) display the zero temperature magnetization 
and susceptibility, respectively for (N) = 1000. Panel (c) 
depicts Curie temperatures versus n c /n;. The filled circles 
correspond to the undamped (I — oo) case and the open cir- 
cles to the damped (I = a) RKKY model. The inset of (c) 
shows VCA and MC T c 's for I = a. Panel (d) is a graph of T c 
versus Xi with n c /rii — 0.03 for / = oo and I = a. For I = oo, 
open circles are Monte Carlo results while the solid line in 
(d) is a theoretical curve, given by 29. 2x 4 J A ; for I — a, closed 
circles are Monte Carlo T c 's, and the solid line is a guide to 
the eye. 

On the other hand, for larger n c /n>i it is also not evi- 
dent in the graph of the magnetization in Fig. 1 (a) where 
the polarization drops to zero (as would happen in the 
thermodynamic limit), and one has instead a long tail 
for n c /rii > 0.1. This slow decay of the magnetization 
for even fairly large systems (e.g. (N) = 1000) discour- 
ages the simple approach of locating the FM/NF phase 
boundary by seeking where m vanishes, and we instead 
seek a sharper signal for the phase transition, magnetic 
percolation. In graph (b), the T = magnetic suscepti- 
bility x = [(S 2 )] - [(|S|)] 2 (with a peak at n c /n t ~ 0.1) 
is shown. Similar singular behavior is also a feature of 
the finite temperature second order FM to NF transition 
(where the temperature is varied with all other parame- 
ters held fixed), and the peak in x(T = 0) suggests that 
our FM to NF phase transition represents a zero temper- 



ature critical point where the critical behavior is driven 
by disorder rather than by thermal fluctuations, a hall- 
mark of a percolation transition. 

Panel (c) of Fig. 1 depicts the Curie temperatures T c 
as a function of n c /rii\ for the undamped case (filled 
symbols) the T c curve attains a maximum value for 
n c /ni ~ 0.03 and then declines, vanishing for n c /n, w 
0.1 in sharp contrast with Curie- Weiss continuum mean 
field theory, which yields the monotonically increasing 

1 /3 

T c oc 7iin c ■ For the damped RKKY model (open sym- 
bols) with I = a, the ferromagnetic transition temper- 
ature peaks for n c /rii ~ 0.12, and is nonzero over a 
broader range than for the pure RKKY model. However, 
as can be seen from the graph, the highest T c for the 
damped model is considerably suppressed relative to the 
peak Curie Temperature for the I — oo case. Panel (d) of 
Fig. 1 also displays T c , but with n c /rii = 0.03 fixed and 
the impurity concentration Xi allowed to vary. The open 
circles are the Monte Carlo results for I = oo, while the 
solid curve is a theoretical curve given by T = 29.2xf^ 3 . 
It can be shown that the good agreement of the x^ 3 law 
with the calculated Curie temperatures implies that even 
for Xi as large as 0.2 the system is in the dilute limit (i.e. 
thermodynamic quantities are insensitive to the details of 
the lattice structure). However, this does not mean that 
the dependence of T c is correctly given by continuum 
MFT, and one actually has T c = x^ 3 'g{n c /rn), where g 
is constant in continuum Curie- Weiss theory, but in our 
case has a nontrivial dependence as highlighted in Fig. 1 
(c) and in the inset, which shows MFT and Monte Carlo 
T c 's on the same graph for I = a. 

By finding the critical n c /rii values for various Mn 
concentrations Xi, we construct a phase diagram for the 
RKKY model for DMS systems, and the result is shown 
in panel (a) of Fig. 2 for the undamped (/ = oo) case. 
From the vertical axis, one sees that a substantial range 
of Xi values is included, certainly encompassing the ex- 
perimentally relevant range. Nonetheless, there is little 
variation of location of the phase boundary which ap- 
pears essentially as a vertical line at the low (but finite) 
n c /rii = 0.1. This is consistent with the notion that the 
dilute limit has been reached even for Mn doping levels 
as high as 20%. 

As ferromagnetism in Gai-^Mn^As is found in exper- 
iment to be robust for n c /rii > 0.1, for a more realis- 
tic treatment we consider a damping scale equal to the 
crystal lattice constant, I = a. Since the factor e~ r l l 
attenuates the more distant AF couplings more sharply 
than the ferromagnetic interactions at closer range, one 
expects RKKY damping to extend the FM/NF phase 
boundary to greater n c /rij values, and the broken line in 
Fig. 2, which displays the FM/NF phase boundary for 
ric/rii > 0.25, is consistent with this intuition. It is im- 
portant to note that although damping expands the do- 
main of the ferromagnetic phase by suppressing AF cou- 
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plings more severely than ferromagnetic interactions, the 
fact that the latter are reduced leaves the ferromagnetic 
state more readily disrupted by thermal fluctuations as 
can be seen in Fig. 1 where T c values for the (I — oo) and 
the (I = a) cases are shown together. Note that an alter- 
nate possibility for an extended FM regime in the phase 
diagram could be the Fermi surface warping [4] due to 
band structure effects in GaMnAs. 
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FIG. 2. Phase diagrams for the RKKY model; the solid 
line corresponds to the undamped, full RKKY model while 
the broken line is for the damped case where I = a. 

For Gai-zMn^As, a salient early experimental finding 
is that T c is maximized for Xi = x° pt ~ 0.07 (though this 
is dependent on the details of sample preparation) with 
the Curie temperature decreasing for doping levels above 
or below this optimal concentration. We consider strong 
local superexchange AF couplings (i.e. between near- 
est neighbors on the fee lattice) as a major contribution 
to the weakening of the ferromagnetic phase for higher 
impurity concentrations, where neighboring pairs of Mn 
impurities are more common. These adjacent moments 
in the fee lattice would presumably form spin singlets and 
hence not contribute to the ferromagnetic phase. Since 
the value of the local Mn-Mn coupling J AF is not pre- 
cisely known, we examine a representative set of values, 
and we work in terms of the rescaled superexchange cou- 
pling j AF = J AF /j£ax KY , where J^ky = 4 ^ Jq ig the 
maximum possible ferromagnetic coupling possible (for 
the I = oo case) for two nearest neighbor impurities on 
the fee lattice. In Fig. 3 (a), we show a graph of T c 
versus x for a large AF superexchange (j AF = —10). 
with n c /rii = 0.05 and I = a held fixed; T c peaks for 
x° pt — .06, in reasonable accord with experiment. In 
general for any finite value of J AF , larger Xi values will 
strongly suppress the FM phase due to the direct AF 
coupling between the Mn moments. 

Using the technique employed for J AF = 0, we obtain 
the T = phase diagram for both the damped and un- 
damped cases with results displayed in panels (b) (I = oo) 
and (c) (I = a) of Fig. 3 for several j AF values. In both 
the damped and undamped cases, it is evident that a 
strong to moderate j AF sharply restricts the domain of 
the ferromagnetic phase; in (b) and (c), the ferromag- 



netic region quickly narrows as the impurity concentra- 
tion rises beyond a few percent, effectively cutting off 
ferromagnetism for Xi > 0.1. Only in the dilute limit, 
where neighboring impurity pairs are less abundant do 
the the FM/NF phase boundaries in (b) and (c) revert 
to their position in the j AF = case. For the case of a 
very strong local AF interaction (with j AF = —10), we 
have determined the values of Xi and n c /rii which maxi- 
mize T c . The points where T c is optimized are identified 
in the phase diagrams of Fig. 3 (b) and (c) with open 
circles at the left of both panels. 
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FIG. 3. The graph in panel (a) displays T c versus x for 
n c /rii = 0.05, I = a, and j AF = —10. Panels (b) and (c) 
are phase diagrams for the RKKY model with a local AF 
coupling incorporated for I = oo (b) and I — a (c); the open 
circles at the left denote locations in the diagrams where T c 
is optimized. 

Finally, we comment on the formation of magnetic 
clusters above T c , for the strongly damped case where 
the length scale I of J(r) is much smaller than the 
typical separation l s between spins. It has been sug- 
gested [10] that in this limit, one can define a cluster- 
ing temperature T* above T c where substantially sized 



4 



magnetic clusters form. In fact, one can argue that 
£ oc Z s (Z s /Z) 2 / 3 [ln(T/T c )]- 2 / 3 , indicating that the clus- 
ter size £ is only weakly dependent on T/T c , with much 
greater sensitivity to l s /l. One can also seek a T* be- 
low which typical magnetic clusters are at least £ = al s 
in size, where a is a dimensionless factor, and one finds 
T* = exp[ci! 3 / 2 (Z s //)]T c ; it is clear that T* rapidly be- 
comes large relative to T c in the I -C l s limit. 

In conclusion, we have worked out the T = phase 
diagram of the the full RKKY model in dilute limit (i.e. 
in the DMS context) at T = 0, finding that a ferromag- 
netic ground state is indeed supported, albeit over a very 
small region of the phase diagram, while continuum mean 
field theory erroneously assumes ferromagnetic order as 
the stable T = phase leading to T C MFT oc runj^ , quite 
distinct from our non-monotonic T c . We have found that 
the T = phase diagram of the RKKY model consists 
of ferromagnetic and non-ferromagnetic phases separated 
by a line of T = critical points in which the NF phase 
gives way to ferromagnetic order via the percolation of 
magnetic clusters. We note that a zero temperature tran- 
sition of a NF to a FM phase signaled by magnetic perco- 
lation is not special to the RKKY model in the DMS con- 
text, but is of broad relevance to magnetic systems where 
there is strong disorder and circumstances which allow 
for competing antifcrromagnctic interactions of a tun- 
able strength. We have found that introducing a cutoff in 
the range of the RKKY interaction (arising, for example, 
from disorder or carrier localization effects) and including 
a local AF superexchange term leads to a T = phase 
diagram and T c behavior in reasonable agreement with 
experiment. The FM phase of the RKKY DMS model is 
fragile, existing over a restricted parameter space (I = oo) 
with moderate T c 's or on a broader domain (finite I) of 
parameter space with reduced T c 's. 
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